Vacuum polarization in the spacetime of charged 
nonUnear black hole 



Waldemar Berej* and Jerzy Matyjasek'^ 
Institute of Physics, Maria Curie- Sklodow ska University, 
pi. Marii Curie- Sklodow skiej 1, 
20-031 Lublin, Poland 

Abstract 

Building on general formulas obtained from the approximate renormalized 
effective action, the approximate stress-energy tensor of the quantized mas- 
sive scalar field with arbitrary curvature coupling in the spacetime of charged 
black hole being a solution of coupled equations of nonlinear electrodynam- 
ics and general relativity is constructed and analysed. It is shown that in 
a few limiting cases, the analytical expressions relating obtained tensor to 
the general renormalized stress-energy tensor evaluated in the geometry of 
the Reissner-Nordstrom black hole could be derived. A detailed numerical 
analysis with special emphasis put on the minimal coupling is presented and 
the results are compared with those obtained earlier for the conformally cou- 
pled field. Some novel features of the renormalized stress-energy tensor are 
discussed. 

PACS numbers: 04.70.Dy, 04.62+v 



I. INTRODUCTION 

One of the most intriguing open questions in modern theoretical physics is the issue of a 
final stage of the black hole evaporation. Although a definite answer would be possible only 
with the full machinery of quantum theory of gravity, some important preliminary results 
may be, in principle, obtained within the framework of semiclassical theory. However, as the 
semiclassical approach could not be used to describe the evolution of the system in the Planck 
regime, it can, at best, tell us about a tendency rather than the limit itself. Unfortunately, 
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even this simplified programme is hard to execute as the semiclassical Einstein field equations 
comprise a rather complicated set of nonlinear partial differential equations, and, moreover, 
the source term — the renormalized stress-energy tensor — should be known for a wide 
class of nonstatic metrics. It is natural, therefore, that in order to make the back reaction 
problem tractable, one should refer to some approximations. 

It seems that the most promising approach consists in constructing the approximate 
stress-energy tensor in the geometry of a static black hole and subsequent computation 
of the semiclassical corrections to the classical metric. Although evaluated in the static 
background, such corrections, as has been pointed out in Ref. [I]], are relevant because 
they give direct information about the influence of the quantum effects on the temperature 
of a black hole. Thus far this programme, initiated in Ref. P|, has been carried out for 
massless flelds in Schwarzschild spacetime and for massive scalar flelds with arbitrary 
curvature coupling in Reissner-Nordstrom (RN) geometry where, among other things, 
quantum corrections to the geometry, entropy, and trace anomaly were computed. The most 
important ingredient of the approach is therefore the renormalized stress-energy tensor of 
the quantized fleld propagating in a spacetime of the static black hole constructed in a 
physically interesting state. However, direct evaluation of this object is rather complicated. 

It is believed that the physical content of the theory of quantum flelds in curved spacetime 
is contained in the (renormalized) effective action, Wr, a useful quantity allowing evaluation 
of the stress-energy tensor by means of the standard formula 

— —W = (T"^) il) 

1 /o r ''"ren I ren- \^ j 

g ' dgab 

Unfortunately, the effective action is a nonlocal functional of the metric and its exact form is 
unknown. In the attempts to construct the renormalized stress-energy tensor one is forced, 
therefore, to employ numerical methods or to accept some approximations. 

For the quantized massive flelds in the large mass limit, i.e. when the Compton length is 
much smaller than the characteristic radius of a curvature, the nonlocal contribution to the 
effective action can be neglected, and the series expansion in of the renormalized effec- 
tive action, Wr, may be constructed with the aid of the DeWitt-Schwinger method P|-pT^. 
Thus constructed renormalized effective action has the form 

Wren = [ ^^^^^ E K (x, x')] , (2) 

n=3 ^ ^ 

where [a„(x,a;')] is a coincidence limit of the n-th Hadamard-Minakshisundaram-DeWitt- 
Seely (HMDS) coefficient. The coefficients are local, geometrical objects constructed from 
curvature tensor and its covariant derivatives of rapidly growing complexity. So far, only 
coefficients for n < 4 are known, but it seems the effective action constructed from higher 
order terms would be practically intractable, especially in the attempts to calculate the 
stress-energy tensor. 

The first nonvanishing term of the effective action constructed for the massive scalar field 
with arbitrary curvature coupling satisfying 

(-□ + + m^) (j) = 0, (3) 
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where ^ is the couphng constant and m is the mass of the field, form the coincidence hmit 
of the HMDS coefficient [03] , is given by |12-14|: 
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where 77 = ^ — 1/6 and ctj are numerical coefficients that stand in front of geometric terms 
in Wren ■ Differentiating functionally Wrln with respect to a metric tensor one obtains 
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where each T*^*)"^ is rather complicated expression constructed from the curvature tensor, 
its covariant derivatives, and contractions. Such calculations have been undertaken in Refs. 
T5| , p!6|| , where generic expressions for the first nonvanishing order of the renormalized stress- 
energy tensor were obtained. They generalize earlier results of Frolov and Zel'nikov for 
vacuum type-D geometries [|lT|. One can easily extend this result to the case of spinor and 
vector field as the analogous expressions differ only by numerical coefficients a^. 

It should be emphasised, however, that the above assumptions place severe limitations on 
the domain of validity of the obtained approximation. Especially, it would be meaningless, 
at least in this formulation, to consider the massless limit of the approach. Therefore, the 
result, which consists of approximately one hundred local geometrical terms, may be used 
in any spacetime provided the temporal changes of the background geometry are slow and 
the mass of the field is sufficiently large. Because of the complexity of the thus obtained 
tensor it will be not presented here, and for its full form and technical details the reader is 
referred to |T6| . 

For quantized massive scalar fields with arbitrary curvature coupling in static and 
spherically-symmetric geometries there exists a different method invented by Anderson, 
Hiscock, and Samuel [jl^. Their calculations were based on the WKB approximation of 
the solutions of the radial equation and summation of thus obtained mode functions. Both 
methods are equivalent: to obtain the term proportional to one has to use 6-th order 
WKB, and explicit calculations carried out for the Reissner- Nordstrom spacetime as well as 
in wormhole geometries considered in Ref. yielded identical results. Moreover, detailed 
numerical analyses carried out in Ref. and briefly reported in |I| show that for mM > 2 
(M being the black hole mass), the accuracy of the Schwinger-DeWitt approximation in the 
Reissner-Nordstrom geometry is quite good (1% or better). 
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Geometry of the Reissner-Nordstrom black hole is singular as r — and there is a natural 
desire to construct its regular generalizations. It is expected that a good candidate for the 
source term of the Einstein equations is the (classical) stress-energy tensor of nonlinear 
electrodynamics. Moreover, renewal of interests in nonlinear electrodynamics that is seen 
recently is motivated by the observation that such theories arose as limiting cases of certain 
formulations of string theory. Unfortunately, the no go theorem proved in Ref. [0 clearly 
shows that for Lagrangians C{F), [F = FabF°'^), with the Maxwell weak-field limit there are 
no spherically-symmetric static black hole solutions with a regular center. 

Recently, employing the Schwinger-DeWitt approximation we have constructed the 
renormalized stress-energy tensor of the quantized conformally coupled massive scalar field 
in the spacetime of the electrically charged regular black hole being an exact solution of the 
equations of the nonlinear electrodynamics and the Einstein field equations. Such a solu- 
tion has been proposed by Ayon-Beato and Garcia (ABG) in Ref. pO|]. It should be noted, 
however, the ABG line element is not a solution of the standard nonlinear electrodynamics 
and the effective geometry (i.e. the geometry seen by photons of the nonlinear theory) is 
singular. Fortunately, the ABG solution may be reinterpreted as describing a magnetically 
charged regular solution of the coupled equations of standard nonlinear electrodynamics and 
gravitation with much more regular behaviour of the effective geometry |2T|]. Moreover, it 
has been shown recently that it is possible to combine electric and magnetic line elements 



to obtain a regular electric solution with a magnetic core |22 



For small and intermediate charges as well as at large distances, the geometry of the ABG 
black hole resembles that of RN; noticeable differences appear near the extremality limits. 
It is interesting therefore to analyse how the similarities in the metric structure of the ABG 
and RN black holes are reflected in the structure of the stress-energy tensors. As for the 
conformally coupled massive scalar field the linear, quadratic, and cubic terms in 77 in the 
stress-energy tensor (|^) are absent, it is anticipated that such similarities do occur. Explicit 
calculations carried out in [|TB1 confirmed this hypothesis and showed that for small values 



of q the appropriate tensors are practically indistinguishable, and, as expected, important 
differences appear near and at the extremality limit. 

In this paper we shall extend the results of Ref. |[T^ and investigate a much more com- 
plicated case of arbitrary t] with special emphasis put on the minimal coupling. Although 
complexity of the approximate stress-energy tensor constructed for the ABG line element 
practically prevents its direct examination, it is possible to extract interesting information 
expanding (T^) into a power series and retaining a few leading terms. We shall show that 
such analyses could be carried out for small charges, large distances, and in the vicinity of 
the event horizon of the extremal ABG black hole. Moreover, on general grounds one may 
easily estimate the role played by T^^^""^ and T(3)afe_ gam. a deeper understanding of the 
problem, however, we employ numerical calculations. 

The paper is organized as follows. In Sec. II the essentials of the ABG black hole 
geometry that are neccessary in further development are briefly described. In Sec. Ill 
certain features of the approximate stress-energy tensor of the massive scalar field in the 
ABG geometry are discussed and compared with the appropriate tensors constructed in the 
RN geometry. The results of our numerical analyses are presented in Sec. IV in which we 
discuss behaviour of the components of (T^) in some detail and present it graphically. 
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II. GEOMETRY 



The general Reissner-Nordstrom solution of the Einstein-Maxwell equations describing 
a static and spherically symmetric black hole of mass M, electric charge e, and magnetic 
monopole charge has a remarkably simple form 

ds^ = -A (r) dt'^ + A-^ (r) dr'^ + {dO'^ + sin^ ed(j)'^) , (6) 

with 

A{r) = l-— + —;^. (7) 

Since the charges enter (|^) as a sum of squares the metric structure remains unchanged 
under the changes of charges leaving = + constant. For Q^/M^ < 1, the equation 
A{r) = has two positive roots 

r± = M ± y^W^Q^, (8) 

and the larger one, r_|_, determines the location of the event horizon, whereas the smaller, r_, 
gives the position of the inner horizon. In the limit = horizons merge at r = M and 
the black hole degenerates to the extremal one. At r = the Reissner-Nordstrom solution 
has an irremovable curvature singularity, which, although hidden to an external observer 
for Q^/M^ < 1, remains a somewhat unwanted feature of the solution. The solution for 
Q'^ > is clearly unphysical. 

At least for purely electrical solutions the nonlinear electrodynamics does not remedy 
the situation. Indeed, consider ralization of the nonlinear electrodynamics with the action 
functional of the form 



S = Y^J d'x^g[R-C{F)], (9) 

where i? is a curvature scalar, F = FabF"-^ , and C {F) is an arbitrary function with Maxwell 
asymptotic in a weak field limit, i. e. C (F) F and j^C (F) — >■ 1 as F ^ 0. 
According to a well known theorem [^,^, there are no regular, static, and spherically 
symmetric solutions of general relativity coupled to nonlinear electrodynamics describing 
a black hole with nonzero electric charge. However, as has been explicitly demonstrated 
recently by Bronnikov, adopted hypotheses leave room for appropriate regular solutions with 
a nonzero magnetic charge. In this regard it is interesting to note that within a different 
formulation of the nonlinear electrodynamics proposed in (a V framework according to 
the nomenclature of Refs. P^ , pT|] ) obtained from the standard one (the formulation) 
by means of a Legendre transformation, Ayon-Beato and Garcia constructed a regular black 
hole solution with a nonzero electric charge and mass [^. Their solution has the simple 
form d^), where 

A (r) = 1 - — ( 1 - tanh ^^^j- J . (10) 
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Bronnikov also demonstrated that any spherically symmetric solution constructed within 
the framework has its counterpart with the same metric tensor constructed within the 
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V framework, and therefore the electric solution @ has the magnetic companion with e 
replaced by the magnetic charge e^- Moreover, Burinskii and Hildebrandt have proposed 
recently a regular hybrid model in which electrically and magnetically charged solutions 
were combined in such a way that the electric field does not extend to the center of the 



black hole \22 



Although the magnetic and electric ABG solutions are precisely of the type (H) with (|T0|), 
the geometries seen by photons of the nonlinear theory are different. It is because photons 
of the nonlinear theory move along null geodesies of the effective metric, and the latter is 
singular for the electric solution P5|-p7[]. It should be noted, however, that otherwise the 
physical geometry is regular and well behaving. 

Since geometries outside the event horizon are described by the same line element and 
since we are going to consider neutral scalar field only, our results will hold for any particular 
realization of the ABG black hole as the only concern here is the metric structure of the 
spacetime. Therefore in what follows we denote both electric and magnetic charge by e. 

Inspection of the metric potentials reveals interesting features of the ABG solution: 
its curvature invariants are finite as r — > and at large r it approaches the Reissner- 
Nordstrom solution. Moreover, for small and intermediate values of the ratio e^/M^, the 
radial coordinate of the event horizon, r_|_, is close to the event horizon of the RN black hole; 
significant differences occur near the extremality limit. 



It has been shown in |16] that the location of the horizons of the ABG black hole may be 



expressed in terms of the Lambert function W p8[. Indeed, making use of the substitution 
r = Mx, = q^M^, it could be demonstrated that the location of the horizons is given by 
the real branches of the Lambert functions: 

" ~AWie, -g74 exp(gV4)) - g2 ' 

where e is for the event horizon, x+ , and —1 for the inner one, x_. The functions W{0, s) 
and W{—l,s) are the only real branches of the Lambert function with the branch point at 
s = — e~^, where e is the base of natural logarithms and 

W{Q,-e-^) = W{-l,-e-^) = -1. (12) 

Consequently, for 

q = (13) 

the horizons merge at 

^ 1 + w' ^^^^ 

where w = W (0,e~"^) = 0.2785. Numerically, one has Xextr = 0.871 for |e| /M = 1.055. 

For small values of q, the ABG line element resembles that of RN. Indeed, expanding 
the function A (x) one has 

Aix) = l-l + ^^+0{q^). (15) 
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In spite of similarities between ABG and RN geometries for q <^ 1, there are substantial 
differences in the extremality limit: the extremal RN line element is described by with 

A(x) = (i-iy, (16) 

whereas near the event horizon of the extremal ABG black hole the function A{x) has the 
following expansion 

A{x) = {x-XextrY + ^2w'^ ~ ^^^'^^'^^^ ^ 0{x - Xextrf ■ (17) 

The form of the above expansion indicates that the proper length between the event horizon 
and any point located at r > Textr is infinite. 



III. THE RENORMALIZED STRESS-ENERGY TENSOR 

Structure of terms in the effective potential (^) and in Eq. (^ indicates that for the 
conformally coupled massive scalar field, i.e. for 77 = 0, there are substantial simplifications 
of the renormalized stress-energy tensor as in this very case the third, fourth, and fifth term 
in (D do not contribute to the final result. Moreover, from (^ it is evident that in i? = 
geometries and for arbitrary curvature coupling the functional derivative of the first and the 
third terms in (^ with respect to the metric tensor vanishes and, therefore, the approximate 
renormalized stress-energy tensor of the massive field in RN spacetime has a general form 

n = Cl + vDl (18) 

On the other hand, however, for the ABG metric one has 



R = — %— tanh ( — 
M^x^ V 2x 



tanh 



2x 



(19) 



which for small q is O (q^) . Detailed calculations carried out for g <^ 1 show that for the 
ABG solution neither Ti^^'' nor Ta^^^ contribute importantly to the result. 

Because of the similarities of the metric structure of the RN and ABG solutions, the 
overall behaviour of the renormalized stress-tensors for conformal coupling should be qual- 
itatively similar and comparable in magnitude at least for small and intermediate values 
of q. However, the differences between the line element ([T6| ) and the expansion (^^), give 
strong evidence for differences in the appropriate components of the approximate stress- 
energy tensors. For the conformally coupled massive scalar field this statement has been 



confirmed by extensive numerical calculations reported in Ref. |]T6[. Now we shall analyse 
the case of an arbitrary coupling with curvature. First, observe that contribution of t^^T^^^"^ 
and 7^3j'(3)a6 ^Yie stress-energy tensor could be made arbitrarily great by a suitable choice 
of the conformal coupling. It should be noted, however, that such great values of rj are 
clearly unphysical and should be rejected. Therefore, as the particular case of the conformal 
coupling has been considered earlier, here we shall confine ourselves mostly to the important 
and physically interesting case of minimal coupling rj = —1/6. 
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In order to construct and D'^ one has to either calculate the curvature tensor and 
its covariant derivatives to required order for the line element (P) with and employ the 
general formulas presented in |]T5[ or to make use of the method proposed by Anderson, 



Hiscock and Samuel [|T^]. Both approaches yield the same result, which reads: 

t _ 810 xY - 855 - 202 x^g^ + 1878 - 1 152 x^g^ _ 2307 x^q'^ + 3084 xg^ - 1248 
* ~ 30240M67r2m2xi2 ' 

(20) 

^, - 1008 x=^ g2 ^ 819 g^ + 2604 xV + 728 xV- 2712 xg^ + 360x^- 792x3 

1)'^ = ; ; ; ; t 1 91 i 

* 720M^7T^mV ' ^ ^ 



^ _ 842 x^g^ + 444 q^ + 162 x^g^ - 462 x^ - 1488 x'^g^ - 1932 xg^ + 315 x^ + 2127 x^g^ 

C+ — . „ , ,^ „ — TT—r^ , 

(22) 



't 30240M6-2^2^i2 



_ 504 xg^ - 208 x^g^ - 588 x^g^ + 336 x^q^ - 117 g^ + 216 x^ - 144 x^ 

720MH^m^x^^ ' ^ ^ 

g _ 2202 x^ - 486 x^g^ - 945 x^ - 3044 x^g^ + 4884 x^g^ - 9909 x^g^ + 10356 xg^ - 3066 g^ 
~ 30240M67r2m2xi2 

(24) 

and 

g _ -1176 x3g2 + 1053 g^ + 3276 x2g2 + 832 x^g^ - 3408 xg^ + 432 x^ - 1008 x^ 

^ ~ 720M67r2m2xi2 ' ^ ^ 

The analogous tensor in the ABG geometry is more complicated, and besides the linear 
terms it contains also terms which are quadratic and cubic in rj. Each component of the 
stress-energy tensor consists of more than three hundred terms and has a general form 

„ „ , ^„ I 1 — tanh — I Qi-j-fc— tanh^ — , (26) 
967r2m2M6 V 2x J ^ ^ x^ 2x ^ ' 



where aijk for a given rj are numerical coefficients (0 < i < 6, 8 < j < 15, < A; < 8), and 
for obvious reasons will not be presented here. 

Since the form of the Eq. (^) differs considerably from the stress-energy tensor given 
by (0) with (p0|-p5|), it could be expected that the runs of both tensors have nothing in 



common. On the other hand similarities in the metric structure of the RN and ABG black 
holes discussed earlier strongly suggest the opposite. Unfortunately, the complexity of the 
renormalized stress-energy tensor of the massive scalar field in the ABG geometry practically 
invalidates direct analytical treatment. One can, however, obtain interesting and important 
information analysing certain limiting cases. Indeed, expanding the stress-energy tensor for 
g -C 1 one has 



8 



^j.yBG ^ J., + ^(1). + (^8) ^ (27) 

where is given by (O-pSf) and 



(28) 



2 2A{i)t_ gMl2049 - 11416 x + 2660 a;2) 2 (617 - 596 x + MOx^) 
12 g6 (459 - 366a; + 70a;2) ry^ 



M6 

2 2A(i)r _ _ (1555 - 1064 X + 140 x^) _ 2 (121 - 140 x + 40 x^) r] 
" ~ 420 M6 xi2 5M6xi2 

12g6 (81 -84x + 20x2) 

+ ]^6^2 ' (29) 

2 2Aa)e _ (5249 - 3528 x + 560x2) ^ 2 (779 _ 720 x + 160 x^) r/ 



(30) 



420 M6 xi2 5 M6 a;i2 

12 g'^ (540 -423 x + 80x2) ^2 



M^x 
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Inspection of expansions (P8H30|), which are vahd for any rj and x, shows that in this order 



the term T^g''^ is absent in the final resuh. It is because Ta^^'' ~ O (g^^) • The second term 

which vanishes in the RN geometry is now O (g^) . 
A similar expansion valid for large values of radial coordinate has the following form 

^j.yBG ^ j.b ^ ^(1)6 _^ ^(2)6 ^ (^-13) ^ (31) 

where is given by (0^, Ai^^*" by (|8|-|g), and 



9W2W g* (-2574 + 1015 g2) Mg^ (-39 + 20g2) « 105 g^ (-27 + 8 g2) 772 
* 210M6xi2 5M6xi2 ^ M6a;i2 

(32) 

2 2 a (2). _ (-738 + 35 g2) 4 g^ (-51 + 28 g2) 5 g^ (-333 + 112 g2) 7^2 

378M6xi2 + 9M6xi2 3M6xi2 ' ^'^'^^ 

^.2 2 A (2)9 gMl2456 - 875g2) 112 g^ (-12 + 5 g2) 5 g^ (-1989 + 560 g2) r/2 
^ 1890M6xi2 9M6xi2 ^ 3M6xi2 

(34) 

Note that in Eqs. (|32|-p^) there is no contribution from the term Ta and, consequently, 
the result which is valid for any combination of couplings and charges is independent of 77^. 

In the vicinity of the event horizon of the extremal black hole the renormalized stress- 
energy tensor has the following expansion 

967r2m2(T,^)^^^ = + t(2)^ (x - x,,,,) + O {x - x^.*.)' , (35) 
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if — ij. 



1 {l + wf{2 + w)(w-lfr]^ 1 {l + wf{4: + w + w^ + 2w^)r] 



f(2)t _ .(2)r 
''t ~ ''r 



4096 


w^M^ 245760 


w6M6 




1 {l + w] 


^(w;=' + 3w2 + 3w + 5) 




(36) 


' 1935360 






3 {l+wY{w 


+ 3)(«;2-2w + l) r^^ 1 


{l + wY {w^ 


- 3) + 3) 



16384 2580480 w^M^ 

1 (l + «;)'(to + 3) (3ti;^-2^; + 7)r7 
983040 ' ^ ^ 

1 {l + wf{w + b){w-lfr]^ , 1 (l + «;)'(-^ + 13 + 3^i;2 + w3)^ 



8192 w6M6 491520 w^Me 

1 (l + u;)^(3w + 17 + 3w2 + w3) 



and 



3870720 w;6M6 



(2)0 _ (2),^ 3_ (1 + wY (-25 u> + 34 - 23 + 9 + 5 1/;^) r^^ 

^ ~ ~ 32768 M6w7 



{31 



16384 M^w^ 

1 [l + wY {-185w -190-lllw'^ + 510^^ -15w^) 7] 
^1966080 ■ ^ 



Since Ta vanishes in the hmit x Xextr there are no terms proportional to rj'^ in (|36|^8D. 

The components of the stress-energy tensor of the massive scalar field are regular func- 
tions of the radial coordinate and are finite on the event horizon. Moreover, it could be 
demonstrated that the difference between radial and time components of {T^)^^^ factor- 
izes: 



2M f , \ 

1 1 - tanh * 



r V 2Mr 



) 



F(r), (40) 



where F (r) is a regular function and, consequently, the stress-energy tensor is finite in a 
freely falling frame. To demonstrate this let us consider a slightly more general line element 

^ _f (^1) ^dx'Y + g{x') {dx'Y + {x'Ydn'. (41) 

For a radial motion the vectors of the frame are the four-velocity Cg = and a unit length 
spacelike vector e" = tt.". Then, using the geodesic equation, one finds 



e^o)=-"=|y,-W(y-l)^,0,0) (42) 



and 
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6(1) — ^ 



±y^._^,0.o), (43) 



where 7 is the energy per unit mass along the geodesic. A simple calculation shows that the 
components 7'(o)(o)) ^(o)(i) and in a freely falling frame (independent of the function 

g{x^y) are 

J(o)(o) = J Ji, (44) 

T(i)(i) = ^'^^y^°^ T,\ (45) 

-'(o)(i) ^ • (46) 

One concludes, therefore, that if all components of and 

/ 



(47) 



are finite on the horizon, the stress-energy tensor in a freely falling frame is finite as well. 
It is expected that constructed formulas satisfactorily approximate the exact stress-energy 
tensor. On the other hand however, to establish the regularity of the exact (T^)^^*^ in a 
freely falling frame one has to explicitly demonstrate that (|T7| ) is satisfied. 



IV. NUMERICAL RESULTS 

Considerations of the previous section are limited to analytically tractable special cases. 
However, to gain insight into the overall behavior of the stress-energy tensor for any combi- 
nation of couplings and charges one has to refer to numerical calculations — our complete 
but rather complicated analytical formulas are, unfortunately, not of much help in this re- 
gard. Below we describe the main features of the constructed tensors and present them 
graphically, fixing our attention on the physically interesting case of minimal coupling. Re- 
lated plots showing the radial dependence of C\ and Tf^^ can be found in Ref. [jT6[, where 
the discussion of the stress-energy tensor of the massive scalar field with conformal coupling 
with curvature has been carried out. 

The case of arbitrary 77 is much more complicated as the tensor (T^)"^-^*^, given by (^, 
is modified by the presence of additional terms. However, numerical analysis performed for 
?7 = — 1/6 reveals that the contribution of ifTf^^ is neglegible for q even as large as 1.0. 
Moreover, a closer examination indicates that as one approaches the extremality limit, the 
magnitude of this very term becomes comparable with rj'^Ta'^^^ and rjTa^^'' only in the closest 
vicinity of the event horizon. Our calculations also suggest that irrespective of the exact 
value of g, each component of Ti^^^ is a monotonic function of the radial coordinate in a large 
neighbourhood of the event horizon, if not in the whole range r > r^. For the extremal case 
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such a behaviour is illustrated in Figs. 1-3, where the time, radial, and angular components 

are displayed. For comparison and completeness we also present the run 
and D'^, where the latter two tensors have been calculated for the RN line element with the 
aid of Eqs. (p0|- |25|) . It should be noted that even in the extremal case, when the differences 
between two considered geometries are most prominent, the tensors Ti^^^ and behave in 
a similar manner though differing noticeably in magnitude. For small values of q, the curves 
constructed for both types of black hole are almost indistinguishable, which could be easily 
established from Eq. (|27|). 



On the other hand however, the components of Ta'^^^ exhibit quite different but very 
regular behavior, which will be described in some detail. For g -C 1 it can be, of course, 
inferred from the approximate formulas (P7|-pil|); greater values of q require numerical exam- 
ination. Specifically, for small charges each component of the considered tensor is negative 
at r = r+, and before approaching zero as r — >■ oo, it changes sign once for the radial com- 
ponent and twice for time and angular components. With the increase of q, all components 
become strongly negative at r+, with their zeros shifted towards larger values of the radial 
coordinate. Near the extremality limit, the tensor Ti^^^ is very sensitive to changes of g, 
which is illustrated in Figs. 4-6. If the value of ratio |e| /M slightly exceeds 1.0, for all 
components there occurs a negative minimum not far away from the event horizon. For the 
extremal ABG black hole Ta vanishes at r_|_, as is clearly seen from Eqs. and 
exhibits oscillatory-like behavior with a rapidly decreased amplitude and increased intervals 
between zeros. The angular component attains a very distinct maximum near the event 
horizon whereas for time and radial component there are minima close to r+. 

The competition of the terms described above results in the overall behaviour of the 
stress-energy tensor of minimally coupled massive fields in the geometry of the ABG black 
hole. Numerical calculations indicate that for small and intermediate values of g, up to about 
0.8, {T^)^^'^ still resembles that evaluated for conformal coupling. For larger values of g, 
however, the change of curvature coupling leads to a considerable modification of the results 
as both magnitudes and radial variations become significantly different. This is illustrated in 
Figs. 7-9, where values of q at and near the extremality limit were chosen. The oscillatory- 
like behaviour of Ti^"*^ for q close to q^xtr is refected by the presence of local extrema visible 
in the stress-energy tensor, and, for the angular component, an inflection point. It is in a 
sharp contrast with the almost monotonic behavior of the stress-energy tensor of conformally 
coupled massive fields depicted in Fig. 7-9 by the dashed lines. Moreover, it should be noted 
that for ?7 = — 1/6, when the extremality limit is approached, there are substantial changes 
of values of {T^)^^'^ at the event horizon as well as in the behavior of the stress-energy 
tensor in a narrow strip near r^xtT- There appear also certain new features for the time and 
radial components which are not present in the case of the conformal coupling. Indeed, for 
0.987 < g < 1.032, the energy density p = —{T^)^^'^ is positive at the event horizon whereas 
for the same values of g, the horizon value of the radial pressure pr = (T^)^^'^ is negative. 
It should be noted in this regard that for 77 = the radial component of the stress-energy 
tensor is always positive there. On the other hand, the angular pressure is positive on the 
event horizon for |e| /M < 0.922 and negative for larger values that is very similar to the 
previously studied case of conformal coupling. 
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V. CONCLUDING REMARKS 



In this work our goal was to construct the renormalized stress-energy tensor of the 
quantized massive field in the spacetime of a nonlinear black hole and to investigate how the 
choice of the curvature coupling affects the results. A regular electrically charged solution of 
this type has been recently proposed by Ayon-Beato and Garcia in V formulation of nonlinear 
electrodynamics and reinterpreted by Bronnikov as a regular magnetically charged solution 
of the standard T formulation. For small and intermediate values of the ratio |e| /M the 
metric structure of the nonlinear solution closely resembles that of Reissner and Nordstrom 
and the similarities in the line elements are refiected in the behavior of the stress-energy 
tensor of the conformally coupled massive scalar fields; notable differences appear near and 
at the extrcmality limit. 

As the general (T^)^^'^ contains terms that are quadratic and cubic in r], the case of 
arbitrary couphng is more complicated. Again, for small and intermediate values of q there 
is a similarity between the (T^) evaluated for the minimally coupled scalar field in the 
AEG geometry and its RN counterpart. Modifications for q slightly exceeding 0.8, which 
are mainly due to the term Ti^"*^, are noticeable although there are still some similarities. 
Comparison of {T^)^^^ for charges close to qextr indicates that the radial dependences for 
both couplings are completely different. Indeed, the behaviour of the stress-energy tensor 
for the minimal coupling is far more complicated as compared with its almost monotonic 
radial dependence for = 0. Moreover, the former is even more sensitive to changes of 
q, especially near the extremality limit. In addition, new features occur as, for example, 
positivity of the energy density at the event horizon for certain values of q. 

It should be emphasized, that being local the Schwinger-DeWitt approximation does not 
describe particle creation, and, therefore, it must not be applied in strong or rapidly varying 
gravitational fields. However, it is expected that for sufficiently massive fields the method 
provides a good approximation of the exact renormalized stress-energy tensor. 

Finally, we remark that although complicated, the derived stress-energy tensor may be 
employed as a source term of the semi-classical Einstein field equations. Indeed, preliminary 
calculations indicate that it is possible to find an analytical perturbative solution to the 
back reaction, and what makes this issue even more interesting and worth further studies is 
the regularity of the geometry of the ABG black hole. We intend to return to this group of 
problems elsewhere. 
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FIGURES 




FIG. 1. The radial dependence of the rescaled time component AT^ (A = 967r^M^m^) (solid 
thick line) and \T^^^^ (solid thin line) contributing to the renormalizcd stress-energy tensor of the 
massive scalar field for the extremal ABG geometry. The dashed lines correspond to the appropriate 
tensors for extremal RN black hole: \Cl (thin curve) and \D\ (thick curve). 



15 




-0.1 



FIG. 2. The radial dependence of the rescaled radial component AtJ^'*^ (A = ddir'^M^m?) 
(solid thick line) and XTr (solid thin line) contributing to the renormalized stress-energy tensor 
of the massive scalar field for the extremal ABG geometry. The dashed lines correspond to the 
appropriate tensors for extremal RN black hole: (thin curve) and XD^. (thick curve). 
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FIG. 3. The radial dependence of the rescaled angular component XT^^^^ (A 



(solid thick line) and XT^^^^ (solid thin line) contributing to the renormalized stress-energy tensor 
of the massive scalar field for the extremal ABG geometry. The dashed lines correspond to the 
appropriate tensors for extremal RN black hole: XCg (thin curve) and XD^ (thick curve). 
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FIG. 4. The radial dependence of the rescaled time component AT^^^^* (A = OGvr^M^m^) con- 
tributing to the renormahzed stress-energy tensor of the arbitrarily coupled massive scalar field in 
the ABG geometry. The thick curve corresponds to the extremal case and the thin curves are for 
a series of q values approaching the extremality limit, g = 0.9,l,1.05. 
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FIG. 5. The radial dependence of the rescaled radial component XTr (A = 967r^M^m^) 
contributing to the renormalized stress-energy tensor of the arbitrarily coupled massive scalar field 
in the ABG geometry. The thick curve corresponds to the extremal case and the thin curves are 
for a series of q values approaching the extremality limit, q = 0.9, 1, 1.05. 
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FIG. 6. The radial dependence of the rescaled angular component XTg (A = 967r^ M^m^) 
contributing to the renormalized stress-energy tensor of the arbitrarily coupled massive scalar field 
in the ABG geometry. The thick curve corresponds to the extremal case and the thin curves are 
for a series of q values approaching the extremality limit, q = 0.9, 1, 1.05. 
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FIG. 7. The radial dependence of the rescaled time component X{Tf)^^'^ (A = 967r^M^m^) 
of the renormahzed stress-energy tensor of the minimally coupled massive scalar field in the ABG 
geometry solid lines) as compared to the case of conformal coupling in this geometry (dashed lines) . 
The thin curves are for q = 1.02 and the thick curves are for the extremal case. 
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FIG. 8. The radial dependence of the rescaled radial component X{T^)^^'^ (A = 967r^M^m^) 
of the renormalized stress-energy tensor of the minimally coupled massive scalar field in the ABG 
geometry (solid lines) as compared to the case of conformal coupling in this geometry (dashed 
lines). The thin curves are for q = 1.02 and the thick curves are for the extremal case. 
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FIG. 9. The radial dependence of the rescaled angular component X{Tq)^^'^ (A = 967r^ M^m^) 
of the renormalized stress-energy tensor of the minimally coupled massive scalar field in the ABG 
geometry (solid lines) as compared to the case of conformal coupling in this geometry (dashed 
lines). The thin curves are for q = 1.02 and the thick curves are for the extremal case. 
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